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magnetic fields from cosmic microwave background 

Maresuke Shiraishi 1 ' 2 '!^ and Toyokazu Sekiguchi 2 '[|J 

Dipartimento di Fisica e Astronomia "G. Galilei", 

Universitd degli Studi di Padova, via Marzolo 8, 1-35131, Padova, Italy 

^Department of Physics and Astrophysics, Nagoya University, Aichi 464-8602, Japan 

(Dated: April 30, 2013) 

Primordial magnetic fields (PMFs) create a large squeezed-type non-Gaussianity in tensor pertur- 
bation, which generates non-Gaussian temperature fluctuations in the cosmic microwave background 
(CMB). We for the first time derive an observational constraint on such tensor non-Gaussianity from 
observed CMB maps. Analyzing temperature maps of the WMAP 7- year data, we find such tensor 
non-Gaussianity is consistent with zero. This gives an upper bound on PMF strength smoothed 
on 1 Mpc as B\ m P c < 3.2nG at 95% CL. We discuss some difficulties in constraining tensor non- 
Gaussianity due to spin and angle dependence of resultant CMB bispectrum. 



I. INTRODUCTION 

Primordial non-Gaussianity is a powerful probe of in- 
flationary models and various aspects of its property, 
e.g., amplitude and scale dependence, have been inves- 
tigated from a diversity of cosmological and astrophysi- 
cal observables. To date, methods to estimate parame- 
ters characterizing non-Gaussianity in primordial pertur- 
bations from the cosmic microwave background (CMB) 
have been extensively investigated by many authors [lj- 
8]. Some specific types of non-Gaussianity have already 
been constrained by observed data, e.g., f$£ = 2.7 ±5.8, 



NL - -42 ±75 and /^ h = -25 ±39 (68% CL) |9j (For 



/■cq 
J] 



scale-dependent non-Gaussianities, see Ref. |io| ) 

These bounds have been estimated under an assump- 
tion that the primordial non-Gaussianity arises from the 
scalar perturbations. On the other hand, there exist var- 
ious models for the early Universe which predict non- 
Gaussianities associated with not only scalar mode but 
also vector and tensor modes [llT - tlq ] . Despite that many 
attempts have been made so far to constrain primor- 
dial non-Gaussianities, those in vector and tensor per- 
turbations predicted by these models have yet to be 
constrained so far. Provided predictions from theoret- 
ical models and precise data from current CMB obser- 
vations, we believe it is timely to investigate constraints 
on non-Gaussianities in perturbations other than scalar 
ones. Among various theoretical models, we in this pa- 
per focus on the electromagnetic field in the early Uni- 
verse as a mechanism to generate vector and tensor non- 
Gaussianities (171, Il8j. 

By cosmological observations of galaxies, cluster of 
galaxies and cosmic rays, the existence of large-scale mag- 
netic field at the present Universe is supported (See e.g. 
Refs. [19|, l20|)- There have been a number of studies in 
which vector fields which exist during inflation are exam- 
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ined as sources for the observed magnetic field |2lTl24| . 
However, so far there exist no viable models for the mag- 
netogenesis via primordial vector fields [251429] . While 
these theoretical considerations strongly restrict model- 
building, phenomenological approaches to constrain pri- 
mordial magnetogenesis are also important. Specifically, 
through impacts on the CMB anisotropy, properties of 
primordial magnetic fields (PMFs) which are assumed to 
be generated from primordial vector fields can be con- 
strained. For example, observational constraints on the 
amplitude of PMFs as well as its scale-dependence can be 
obtained by the CMB power spectra alone (For current 
bounds, see e.g. Refs. |30l - l33l |). 

Assuming that field strength of PMFs has a Gaus- 
sian distribution, its anisotropic stress fluctuation cre- 
ates both the scalar and tensor squeezed-type non- 
Gaussianities due to the quadratic dependence on the 
field strength. This leads to non-Gaussian CMB 
anisotropy and suggests that higher order correlation 
functions or polyspectra of CMB anisotropy beyond the 
power spectrum should also be informative in probing 
PMFs. On the basis of this concept, this paper newly 
explore an observational constraint on the PMF strength 
by evaluating the magnitude of non-Gaussianity in the 
CMB temperature anisotropy. 

In the case of PMF, the tensor non-Gaussianity domi- 
nates over the scalar one [34J and hence the non-Gaussian 
temperature fluctuation mainly holds information of the 
tensor mode. Since CMB fluctuations generated from 
PMFs have unique features and are distinct from those 
from ordinary scalar perturbations in inflationary Uni- 
verse, non-trivial constraint can be expected to be ob- 
tained. In this sense, this work corresponds to a first at- 
tempt to constrain tensor non-Gaussianities from CMB 
data. This is also another motivation of this paper. 

This paper is organized as follows. In the next sec- 
tion, we discuss how to simulate non-Gaussian temper- 
ature maps generated from PMFs. In Sec. IIII1 we esti- 
mate limits on tensor non-Gaussianities from PMFs using 
the observed temperature maps of the WMAP 7-year re- 
sult and translate them into observational bounds on the 



PMF model. The final section is devoted to summary 
and discussion. 



II. NON-GAUSSIANITY IN THE CMB MAP 

GENERATED FROM PRIMORDIAL MAGNETIC 

FIELDS 



First, we shall briefly summarize the mechanism of 
PMFs to generate CMB temperature fluctuations and its 
signatures in the observed CMB maps. After PMFs are 
produced and stretch beyond horizon during inflation, 
the anisotropic stress of PMFs acts as a source term in 
the Einstein equation and supports the growth of cur- 
vature and tensor perturbations even on superhorizon 
scales until neutrino decoupling. However, subsequent 
to neutrino decoupling, finite anisotropic stress fluctu- 
ation in neutrino cancels out the magnetic anisotropic 
stress fluctuation and therefore the enhancement of the 
metric perturbations ceases. The resultant curvature 
and tensor perturbations produce the CMB scalar and 
tensor anisotropies, which are called passive mode fluc- 
tuations [33. Supposing that the PMFs are quantum- 
mechanically created and the probability distribution of 
their field strength obeys pure Gaussian statistics, the 
PMF anisotropic stress fluctuation becomes a highly non- 
Gaussian quantity due to the quadratic dependence on 
the Gaussian PMF strength. This anisotropic stress fluc- 
tuation has the local-type shape and therefore the bis- 
pectrum of resultant metric perturbations is amplified 
in the squeezed limit. Consequently, CMB anisotropies 
generated from PMFs should be non-Gaussian and their 
bispectrum can be nonzero [U GJ, Hj. The CMB 
temperature anisotropy for given direction n is quanti- 
fied via the spherical harmonics expansion as — J^- = 
Efa fl '"i^™(")- ^ the rotational invariance is pre- 
served, the CMB bispectrum is expressed as 



should be given as 



n 
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Note that B^ 2 ^ z is proportional to the magnetic field 
strength to the sixth power. According to Ref. [34]], in 
the case of the PMF, the tensor perturbation dominates 
over total signal of the CMB bispectrum for £ < 500, and 
the contributions of the scalar and vector perturbations 
are small. Therefore, in this paper, we consider that the 
CMB bispectrum in Eq. (JTJ) is generated only from tensor 
non-Gaussianity. 

To extract the information of the PMFs from the ob- 
served CMB maps, we need to simulate a,£ m from the 
theoretical CMB bispectrum. According to Refs. |36l - 
l38l |. given a power spectrum Ci and bispectrum Bi x i 3 i 3 , a 
random realization of CMB temperature anisotropy a,£ m 
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Here, af m is the Gaussian part of the realiza- 
tion whose variance is given by (rin=i a 2 l ?7 

Ce 1 (— l) mi S£ 1 .£ 2 5 mi ,- m2 , and af® denotes the non- 
Gaussian part of the realizations. 

Eq. §5$ indicates that to generate a single realization, 
a direct implementation requires C(£^ ax ) arithmetics, 
where £ max is the maximum multipole. For £ max ~ iOOO, 
required computational time is enormous. In the litera- 
ture [3a - l38| , the reduced form of a^ has been found for 
the case of the standard scalar non-Gaussianity where the 
angle-dependence is removed. In the CMB bispectrum 
from such non-Gaussianity, the dependence on (li,mi), 
(/?2, 1TI2) and (£3, 1713) can be separated from one another. 
Therefore, in this form, the number of summations can 
be reduced. Unlike the scalar perturbation, the tensor 
perturbation is a spin-2 angle-dependent quantity. In the 
tensor case, this property induces extra correlations be- 
tween £\ , £2 and £3 in the CMB bispectrum, which is not 
seen in the case of scalar bispectrum. In other words, we 
do not have a reduced from of Eq. ^ and hence need to 
straightforwardly perform 0(^f nax ) arithmetics in com- 
putation of ttf£- Fortunately, in the case of the PMF 
model, the signal-to-noise ratio of the magnetic tensor 
mode is saturated at £ ~ 100 and the summations up to 
£ = 100 are enough [3J|, which prevents us from taking 
too much computational time. 

In what follows, we obey the conventional parametriza- 
tion for the power spectrum of PMFs as 

(iFOc^-OO) = (2^) 3 ^pyk)^(k + k'), (4) 
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where P l j(k) = S*j — k l kj, ns and B\ m P c are 
divergence-free projection tensor, the PMF spectral in- 
dex and the PMF strength smoothed on IMpc scale, re- 
spectively. Assuming the generation of the PMFs at the 
GUT scale and nearly scale-invariant shape of the PMF 
power spectrum as ns = —2.9, the lcr error for a cosmic 
variance limited survey with £ max — 100 is evaluated as 
SBi Mpc = 3.1 nG from the Fisher matrix analysis 34 1. 

In the next section, we shall constrain the amplitude 
of bispectrum given by A = {B\ m P c/3.1 nG) 6 oc Bi x i 2 i 3 . 
Note that theoretically, A should take positive values. 
Then, the expected lcr error above corresponds to 5 A = 
1. Later we will confirm that the error estimated from ac- 
tual observational data is slightly relaxed compared with 
this value. 



III. OBSERVATIONAL LIMITS 

Given a theoretical template of B^^^ for a specific 
theoretical model, we can in general construct an optimal 
estimator of the amplitude of the bispectrum of primor- 
dial perturbations [4j, |39| . An optimal cubic estimator 
for the amplitude of bispectrum should be given as 
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where B^ 2 ^ 3 is the template bispectrum for the PMF 
model normalized with A = 1, C^ imii £ 2m2 is the total 
(signal+noise) covariance of observed anisotropy ai m and 
aim = J2£' m '^£m,£'m' a im is a filtered map. A/" is the 
normalization which guarantees the unit response of the 
estimator to true A. In addition, 1/A" corresponds to 
the variance of the estimator in the limit of Gaussian 
perturbations, which therefore gives an error of A. 

Summary of our analysis is as follows. We use the 
foreground-cleaned temperature maps from the WMAP 
7- year observation at V and W bands [4CJ, lU| 1 with a res- 
olution A^ido = 512 in the HEALPix pixelization scheme 
[42J. 2 To reduce effects of residual foregrounds, we apply 
the KQ75y7 mask [Hj] , which cuts 28.4 % of the sky. Fil- 
tered map a£ m is computed based on the method of Ref. 
43]. This method also allows to marginalize over con- 
tributions of Galactic foregrounds adopting spatial tem- 
plates of synchrotron, free- free and dust emissions in Ref. 
[41j as well as the monopole and dipole components. The 
CMB signal power spectrum Ci is computed using the 
CAMB code j44|, assuming a concordance flat power-law 
ACDM model with the mean cosmological parameters 
from the WMAP 7- year data alone j45|. In order to de- 
termine the normalization J\f , we first generate random 
realizations of non-Gaussian CMB maps with A = 1 fol- 
lowing the method we described in Sec. [TTJ After adding 
the simulated noise, we then conduct this mock data into 
the same pipeline as the observed data. The number of 
MC samples we adopt in our analysis is 100, which is 
sufficient for converging results. 

Here we present our constraints. First we include mul- 
tipoles up to ^ max = 100. Without template marginaliza- 
tion of Galactic foregrounds, we obtain A — —5.89 ± 1.06 
at ler, and with template marginalization, A = —2.11 ± 
1.06. This shows that effects of residual Galactic fore- 
grounds are significant at these angular scales. On the 
other hand, if we included multipoles up to £ ma x = 400, 
the constraints change to A = —0.91 ± 1.06 regardless 
for template marginalization of Galactic foregrounds. 3 



Accordingly, we can conclude that the WMAP data is 
consistent with null tensor non-Gaussianity induced by 
PMFs at bispectrum level. Taking into account the above 
result and a theoretical prior, A > 0, we find new upper 
limit on the PMF strength, namely, B\ m pc < 3.2nG at 
95% CL. This is consistent with latest bounds from CMB 
power spectrum data in the SPT or PLANCK experiment 

Mm. 



IV. SUMMARY AND DISCUSSION 



The origin of the observed large-scale magnetic fields is 
one of the most important and interesting issues in probe 
of the early Universe and some researchers seek answers 
in the inflationary paradigm. In this paper, we have dis- 
cussed an observational constraint on the seed magnetic 
field stretched by the inflationary expansion from the 
analysis of non-Gaussianities in CMB anisotropy. Signal 
of the Gaussian field strength of the PMF becomes largest 
on large scales via the enhancement of non-Gaussian ten- 
sor perturbation. 

We have analyzed the data of WMAP 7-year temper- 
ature maps and confirmed no evidence of squeezed-type 
tensor non-Gaussianity due to PMFs. Constraining the 
amplitude of tensor non-Gaussianity leads to an upper 
bound on the PMF strength as B x m P c < 3.2nG (95% 
CL). This result is not sensitive to the residual Galactic 
foregrounds. This value may be improved by considering 
impacts of polarizations or adopting more accurate data 
as the PLANCK data (!]. 

Aside from the issue on the primordial magnetic field, 
this is a first challenge to constrain the primordial tensor 
non-Gaussianity from the CMB bispectrum. The tensor 
CMB bispectrum has a spectral shape quite distinct from 
the scalar one and leads to non-trivial constraints which 
have never seen in the scalar case. Unlike the scalar case, 
the CMB bispectrum composed of the tensor mode has 
the correlation between £i, £2 and £3 due to the spin and 
angle dependence associated with the tensor mode. This 
correlation seems to prevent us from reducing the formula 
for the simulated ag_ m (J3]) to the form with lower compu- 
tational cost, which is derived in the scalar case. This 
may be disadvantage for the data analysis with higher 
resolution and some improvement will remain as a future 
issue. 



1 http://lambda.gsfc.nasa.gov 

http://healpix.jpl.nasa.gov 
3 Note that we here adopt the normalization M from multipoles up 



to £ max = 100. This is because according to Ref. [34| |, we expect 
M should be saturated at (. > 100. While this expectation may 
be affected by the inhomogeneity in noise levels and sky cuts, 
the effects would not be significant since M at £ ma x = 100 shows 
good agreement with the Fisher matrix analysis of Ref. |34l |. 
which assumes homogeneity in noise levels. 
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